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I. Introduction 

An orbiting laboratory provides the unusual scientific opportunity 
to carry out experiments in the absence of gravitational effects. Poten- 
tial benefits in the field of aerosol science are the absence of sedi- 
mentation and convection. In an experiment carried out under terrestrial 
conditions, sedimentation is important for larger areosol particles. 
Preventing spurious convective currents is often troublesome in terres- 
trial aerosol experiments. In addition to these potential benefits, 
the absence of a gravitational field may provide theoretical simplifica- 
tion. The feasibility of an orbital aerosol experiment is therefore 
of scientific interest. This preliminary investigation was undertaken 
to examine the feasibility and scientific benefits for a zero-gravity 
aerosol behavior experiment in an orbiting laboratory. 

II . Ob j ectives 

A. Examine the theoretical feasibility for a zero-gravity 
aerosol behavior experiment. 

B. Examine the technological feasibility for a zero-gravity 
aerosol behavior experiment. 

C. Identify potential scientific benefits for the experiment. 

D. Present conclusions. 

III. Theoretical Feasibility 

In the evaluation of experiments concerned with the disappearance 
of particles from a confined aerosol, one must in general consider the 
simultaneous concentration changes due to coagulation, diffusion, and 
sedimentation. The mathematical complexity of this problem has frustrated 
attempts to devise totally satisfactory analytical models, even 



2 


for initially monodisperse particles in the absence of electrostatic 
and convective effects. The conduct of such experiments under the 
essentially zero-gravity conditions of an orbiting space laboratory 
would provide data in the absence of sedimentation, normally an important 
depletion mechanism for particles larger than a few tenths of one 
micrometer in radius. A key issue is whether the absence of sedimen- 
tation results in a significant reduction in the mathematical complexity 
of the problem. 

Considerable success has been achieved in describing the behavior 
of dilute aerosols in which effects due to coagulation are absent. For 
the. one-dimensional problem, the following partial differential equation 
giving the numerical particle concentration n as a function of time £ 
and vertical position £ has been solved analytically by Davies (1) and 
verified experimentally by Richardson and Wooding (2) : 


3n 

9t 




( 1 ) 


The aerosol considered by these investigators was monodisperse and 
confined between perfectly absorbing parallel surfaces of large hori- 
zontal extension such that effects due to the vertical walls were 
negligible. The particle diffusion coefficient ID and settling velocity 
V may be evaluated from theoretical considerations given by Fuchs (3) . 

A major problem associated with achieving a general solution which 
includes coagulation is the evolution of the particle size distribution 
with time. Hidy and Brock (4) reviewed many investigations of coagulation 
and found that the only exact solution is due to Smoluchowski (5j in 
which the coagulation constant is independent of time. In the absence 
of particle depletion due to diffusion and sedimentation, the basic 
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equation of coagulation is 



( 2 ) 


where IC is 8 ttRD for spherical particles of radius R. Integration of 
equation (2) gives 


i _ i 

n n 

0 


Kt 


■ or 


n _ 1 

n ~ 1 + n Kt 

o 0 


C3) 


where n^ is the initial particle concentration at t = 0. Numerous - 
experimental studies reviewed. by Fuchs, (3) establish the linear relation- 
ship between n”^ and, time although values for j( tend to be somewhat 
larger than computed values. Hidy and Brock (4) suggested that the 
lack of quantitative agreement is attributable to electrical effects, 
spurious air currents, increasing polydispersion, and the need for a . 
slip correction. The effect of increasing polydispersion was investigated 
by Hidy C6) who numerically solved a set of simultaneous nonlinear 
differential equations for a discrete particle size distribution. The 
numerical results give values for n(t3 which closely match those given 
by equation (3) , although changes in the relative concentrations in 
individual particle size categories differ from those for constant 
particle collision parameter. Various authors have suggested that the 
increase in mean particle size to decrease IC is offset by the increasing 
polydispersion with its opposite effect. The linear relationship 
between n ^ and time is used in the present work to develop a macro- 
scopic model for the combined effects of diffusion and coagulation. 

The equation for the general case which includes coagulation gives 


n(x,y,z,t) by 
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II - D - V II - Kn^ (4) 

where v, and K are assiimed independent of time and position. The 
assumptions of perfectly absorbing walls and uniform initial particle 
concentration n^ give n(x,y,z,t) = 0 at the walls and nCx,y,z,0) = n^ in 
the chamber. Equation (4) has been solved analytically by Wilhelm (7) 
for the mathematically similar problem of contained plasma particles 
undergoing simultaneous recombination, diffusion, and convection. 
Wilhelm's approach, which involves a transformation to a differential 
equation in which the nonlinear term becomes a small perturbation, was 
shown to apply to plasma particles confined by nonreflecting walls. 
Application to the coagulating aerosol undergoing simultaneous 
depletion by diffusive deposition was made by Benedict (8) . The 
details of this treatment are given in Appendix I. 

The theoretical feasibility for a zero-gravity aerosol study 
was examined by carrying out simulated experiments with models developed 
by Benedict (8) . The purpose of the computations was to determine 
whether the requirements for a zero-gravity aerosol study are, at 
least in principle, compatible with the time and space limitations 
for an orbital experiment. Input to the models consists of the 
particle properties and initial concentration, gas properties, and 
chamber geometry. Output consists of the particle concentration 
as a function of time and location in the chamber. Numerical 
results are given in Appendix I . 

Two conclusions arise from the theoretical feasibility study. 
Firstly, the time and physical space limitations for an orbital 
experiment are not prohibitive in terms of obtaining kinetic data. 
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For example, experiments of a duration of 1-2 hours in a cylindrical 
chamber (r = 25 cm, h = 50 cm) can produce meaningful data on the 
history of the aerosol confined in the absence of gravitational 
effects. Secondly, the model is internally consistent and produces 
physically reasonable results. 

While the solution procedure resulted in substantial simplification 
of the problem, the reduction in complexity is only partly attributable 
to the absence of the gravitational term. The nonlinear coagulation 
term necessitates a transformation of variables, whether or not the 
gravitational term is present. 

A final point concerns possible limitations of the models 
employed in the theoretical analysis. Both and 1( have been used 
as ensemble parameters. Th,e major limitation of this approach 
is that the models do not provide direct information on the evolution 
of the particle size distribution. It should also be emphasized that 
while equation (1) has been verified experimentally, equation (4) 
must still be regarded as a postulate. The experimental conditions 
under which it is permissable to treat both ^ and 1C as ensemble 
parameters in equation (4) are considered in Appendix I. 

iV. Technological Feasibility 

The technological feasibility was examined by addressing the 
following issues: 

A. Experiment Definition 

B. Requirements for Aerosol Generation 

C. Requirements for Measuring Particle Concentrations 

D. Requirements for Data Analysis 

The major findings are summarized in a previous report (9). Two 
areas of concern have been identified. The first is the current absence 
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of a single, totally satisfactory experimental technique for 

determining aerosol particle concentrations over the particle size 
-7 “3 

range of 10 cm to 10 cm. Either complementary experimental 
techniques would be required or, alternatively, the size range 
accessible must be narrowed to be compatible with a single experimental 
technique. In view of the potential problems associated with 
calibration of complementary techniques under orbital conditions, 
the latter alternative seems more practical. 

The second area of concern is that of determining the particle 
concentration to the accuracy required. For example, in the absence 
of electrostatic and convective effects, calculated and measured 
values of the coagulation constant K differ by perhaps 10%. Clearly, 
the experimental technique selected must be capable of resolving 
differences in particle concentrations smaller than 10% in order 
to make meaningful comparisons between experimental and theoretical 
data. While considerable progress has been made in recent years 
in refining experimental methodologies in aerosol science, the required 
accuracy may not be available for the ranges of particle sizes and 
concentrations of interest. However, recent developments with 
electrostatic classifiers and light-scattering techniques are particularly 
encouraging. The required instrumental capabilities may not 
be far away, but the situation seems borderline at this time. 

V. Scientific Benefits 

Potential scientific benefits of a zero-gravity aerosol study 
include validation of theoretical models for aerosol kinetics and 
measurement of and IC in the absence of convective effects. However, 
because of existing gaps in both the theoretical and experimental 
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aspects, such an experiment would probably be premature at this time. 

In order for the potential benefits to be fully realized, more detailed 
theoretical models for the combined effects of coagulation and diffusion 
are needed. Moreover, additional refinements in certain experimental 
methodologies would be helpful to assure accurate measurements over 
the ranges of particle sizes and concentrations of interest. 

VI . Conclusions 

The theoretical feasibility for a zero-gravity aerosol study 
has been examined. The mathematical complexity of the problem is 
discouraging if one attempts to retain a detailed picture of the 
combined effects of coagulation and diffusion upon the evolving 
particle size distribution. However, an analytical solution is 
possible if one considers only the particle concentration n 
and treats and 1C as ensemble parameters. Experimental studies 
support use of ^ as an ensemble parameter in many cases. Treating 
D_ as an ensemble parameter imposes some limitations, however. The 
results of the macroscopic treatment show that an aerosol decay 
experiment is feasible in a compact chamber for a time duration of 
the order of hours. It is concluded that the limitations of physical 
space and time for an orbital experiment are not prohibitive in terms 
of conducting an aerosol experiment. 

Because of the present mathematical difficulties associated 
with treating the combined effects of coagulation and diffusion, 
reservations are expressed about the scientific urgency for a zero- 
gravity aerosol study at this time. The experiment would also 
appear to stretch existing capabilities for characterization of 
aerosol particles. Nevertheless, the need for reliable and accurate 
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aerosol behavior data in the absence of convective effects is 
recognized. Periodic re-examination of the need for a zero-gravity 
aerosol study is therefore recommended. When such a study is 
planned, it will be important to include a complementary terrestrial 
investigation. The difficulties associated with carrying out an 
experimental aerosol study under favorable terrestrial laboratory 
conditions can be formidable. The success of an orbital experiment 
will be highly dependent upon managing these difficulties and 
anticipating additional problems posed by the orbital situation. 
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ABSTKACT OF THESIS 

THEOBETICAL BEHAVIOR OF A CONFINED AEROSOL 

Coagulation, sedimentation and diffusive deposition are the primary 
removal mechanisms for an aerosol confined in a chamber on earth. The 
equation describing the depletion rate due to these mechanisms is a 
second order nonlinear partial differential equation. For an aerosol 
in a zero-gravity environment the sedimentation term drops out , but 
this does not change the basic nature of the equation. An analytical 
solution to the resulting equation is presented and particle concen- 
trations are computed as a function of time and location in a cylin- 
drical chamber . 

The equation is also solved for an aerosol under the Influence of 
a gravitational field. There are some difficulties with this solution 
w':;are the removal mechanisms are operating at similar rates. These are 
overcome by modelling the decay process as if diffusion were not present. 
Results of the two models indicate that sedimentation is the most impor- 
tant of the removal mechanisms. Coagulation is next in importance and 
diffusion is negligible except within 1 centimeter of the chamber wall. 

Bruce John Benedict 

Mechanical Engineering Department 

Colorado State University 

Fort Collins , Colorado 

Fall, 1977 
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CHAPTER r 
EQUATION 
Introduction 

The principal removal mechanisms for an aerosol confined in a 
chamber which is located in a gravitational field are sedimentation, 
coagulation, and diffusive deposition at the wall. In most situations 
of interest sedimentation is the most important. An experiment carried 
out in the zero-gravity environment of an orbiting space craft would 
provide a unique opportunity to study aerosols. By placing an aerosol 
in such an environment it is possible to obtain measurements of the 
coagulation coefficient and diffusion coefficient of the particles. 

This study provides a model with which such an experiment can be de- 
signed. A model which describes the depletion of an aerosol stored in a 
gravitational field is also presented. 

For the purposes of this modelling effort the following assumptions 
were made: 1. The particles are unit density spheres. 2. The system 

is monodlsperse. 3. Initially, the particles are uniformly distributed 
within the chamber. 4. The coagulation coefficient is constant with 
respect to time. The particle formed when two particles stick together 
should have a coagulation coefficient which is larger than that of a 
single particle because the combined particle is larger. However, Fuchs 
[1] states that this effect is almost balanced by the decreased diffu- 
sivity of the larger particle. 5. All collisions, particle-particle 
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and particle-wall, have a sticking coefficient of unity. Fuchs [1] 
indicates that as long as the only driving forces in the system are 
thermal, the sticking coefficient is nearly one. 

Equation 

Richardson and Wooding [ 2 ] present an equation describing the 
depletion of a monodisperse aerosol which is confined in a chamber in a 
gravitational field. 

Il . DV“n - v|f - Kn^ (1) 

where n is the numerical concentration of particles, t is time, D is the 
particle diffusion coefficient, v is the terminal settling velocity of 
the particle, z is the vertical coordinate of the chamber, and K is the 
particle coagulation coefficient- The boundary conditions are implicit 
in the assumptions. Because the sticking coefficient for particles 
colliding with walls is unity, a particle which strikes a wall sticks to 
it and is thus removed from the bulk of the chamber. Thus we have 
n(r,z,t) = 0 at the walls. The initial concentration was assumed to be 
uniform throughout the chamber, thus the initial condition is n(r,z,0) = 

n . 

o 

Coefficients 

Strauss [3] gives a procedure for calculating the terminal settling 
velocity for aerosol particles in terms of the particle Re 3 molds number, 
R^, and the particle drag coefficient, C^- The procedure is to first 
calculate the product, using properties of the particles and the 

fluid in which they are dispersed. 
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C,R^ 
d e 


32p(p^ - p)R^g 


( 2 ) 


where p is the fluid density, p^ Is the particle density, g Is the 
acceleration of gravity, R Is the particle radius and y Is the viscosity 
of aly. Davies [4] gives a series of empirical formulas for the Reynolds 
number In terms of the product Once the Re3molds number Is deter- 

mined the product can be used to find a value for the drag coefficient. 
These values are then used In Stokes law to find a velocity, V, 


16Rpg(p - p) 

V E (3) 

The terminal settling velocity is determined from V by applying the 
Cunningham slip correction, 


C = 1 + |(1-257 + 0.400exp(-l.lo|)) 


( 4 ) 


where C is the Cunningham slip correction, 6 is the mean free path of 
the fluid molecules. Now It Is possible to determine a value for the 
terminal settling velocity, v. 


V » CV (5) 

The diffusion coefficient Is given in Fuchs [1] 

D * kTB C6) 



20 


where k is Boltzmann’s constant, T is absolute temperature, and B is the 
particle mobility also given by Fuchs IlJ. 


6irR,y 


( 7 ) 


The constants A, Q, and b are given in Millikan [5] for oil drops in air 
at 23Q and 1 atmosphere. They are A - 0,864, Q = 0.290, and b = 1.25. 
The coagulation coefficient is in Fuchs [1] as follows: 


K = 8ttRD 


( 8 ) 


Solution 

The nondimens ional variables, r' = r/R*, z' = z/R*, n' = n/n , and 

o 

t* = tD/R*^ are introduced. R* is the radius of the chamber. 
Substituting these into equation (1) results in 


9n’ 

at’ 


V^n' 


'3n ’ 
^3z’ 




12 


(9) 


where 


_ R*v 
D 


Y = 


R*^Kn 
— — o 


( 10 ) 


The boundary conditions become n^(r',z',t') = 0 at the walls and 
the initial condition is n^(r', 2 ^, 0 ) = 1. 

The nonlinear coagulation term in eqjiation (9) presents the 
greatest difficulty in obtaining a solution thus it is desirable to 
remove it using a transformation. This is done by using the solution of 
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the Smoluchowskl coagulation equation as was done by Wilhelm [6]. The 
cransformatlon used was 


n'(r',z',t’) 


uCr’.z’.t*) 

1 + Yt'u(r',z’,t’) 


( 11 ) 


This results in the equation 


3u 

at' 


= V^u - & 


3u 

3z' 


e 


( 12 ) 


where 


e = 


2yt * r 


■)' 




(13) 


The boundary conditions* remain the same for the transformed variable, 
u, as for the variable, n. That is u(r',z',t') = 0 at the walls and 
u(r' ,z' ,0) = 1. 

The times of 1 and 2 hours, which are used in the calculations, 

are very short as compared to the characteristic time for diffusion. 

—8 

This is bom out by the values of t', which range from 7.0 X 10 for 

-3 

large particles to 1.7 X 10 for small particles. The values of y 
relatively large, ranging from 15 for small particles to 1.6 X 10^ for 
large particles. However, the values are such that for a given particle 
the product yt' is of order 1 or smaller. It is further anticipated 
that except near the walls the spatial derivatives will be small. Thus 
the value for e will be small. 

Equation (12) can now be solved by a successive approximation 
technique. To do this £ is set equal to 0 and equation (12) is solved. 
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This Is the zeroth approximation. Using the solution for the zeroth 
approximation a value can be calculated for c. This value is then put 
back into equation (12) as a constant and the equation Is solved, 
resulting in the first approximation. This process is repeated using 
the value of e calculated from the previous approximation until the 
desired accuracy is obtained. 
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CHAPTER II 

ZERO-GRAVITY SOLUTION 
Equation 

The simplest case for obtaining a solution is the zero-gravity 
situation. This is because with zero-gravity the terminal settling 
velocity of the particle is zero, thus 3=0 and equation (12) reduces 
to 


9u 

3t’ 


= V^u - e 


(14) 


Setting e = 0 for the zeroth approximation results in the diffusion 
equation which for a cylindrical chamber is 


9u _ 1 9u 9^u I 9^u 

3t' " r’ 8r' 97^ 


(15) 


Solution 

The solution for equation (15) can be found using the technique of 
Separation of Variables as given by Wylie [7]. For this it is assumed 
that the solution is the product of three functions, a function of r’ 
only, a function of z’ only, and a function of t* only. Thus 
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u(r',z',t*) = U^(r')U^(z’)U^(t’) 


(16) 


Substituting this into equation (15) gives 


au^ 

at’ 

u 


au 

^ u 

r 


a^u 

r 

aF^ 

u 


a^u 


+ 3z 


z 

TTT 


U 


= V 


(17) 


One side of this equation is a function of t’ only, while the other is 
a function of r' and z\ The only way this can be true is if both sides 
are equal to a constant, V- Using this fact it is possible to solve 
for U^. That solution is 


= exp (vt * ) 


(18) 


This equation implies that V < 0 because the solution for u(r\z^,t') 
cannot increase exponentially with time. Now it is possible to take 
the part of equation (17) which is a function of r’ and and rearrange 
it so that there is an equation which is a function of r' equal to a 
function of z’ which again must be equal to a constant. 


3U 

I 

3r’ 


9^U 


+ 3r 


r 

TT 




3^U 

z 

u 


= a 


(19) 


Taking the z* equation it is possible to solve for U 


U = sin(y^z^) 
z 


( 20 ) 
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From the boundary conditions on the ends of the chamber, a can be de- 
termined . 


a = 


'' H ^ 


( 21 ) 


where H is the height of the chamber. 

The solution of the remaining equation is U^. 


U 

r 


^ J 

o 


(A^r’) 


( 22 ) 


where the are the zeros of >Jq(A^) = 0. Combining all of these solu- 
tion results in the complete solution results in 


u(r',z’,t') = ^|^A^^J^(X.r')sin(^^’)exp(-(X2 + 

(23) 


Using the initial condition, can be found. 


/ V^r * J (A .r ’ )sin("^^^^ * )dz ’dr ' 

^ Q 0 ^ i M 

/ V^r ’ (X . r ’ ) sin^ ( ^~ z ’ )dz Mr ' 
0 0 o 1 H 

^ 4(1 - (-1)^) 
jlTX^J^(X^) 


where L * H/R*. 

Inner-Outer Expansion 

Results of the solution given by equation (23) are that u(r’,z’,t*) 
is almost a constant for most of the chamber, indicating that coagula- 
tion is much more important as a removal mechanism for the bulk of the 
aerosol than is diffusive deposition. The fact that the higher order 
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approximations are very difficult to solve indicates that the solution 
might be found using an inner-outer expansion technique. The inner- 
outer expansion technique is used on problems which have several 
regions in which the processes occurring are of different nature. The 
classical example of such a problem is a boundary layer problem such 
as illustrated here. In this technique a solution is found describing 
the dominant process in each region. The boundary conditions for the 
solution in one region are then modified by using the boundary condi^ 
tions which were shown to exist in the adjacent region. These modifi- 
cations are repeated until the boundary conditions for adjacent regions 
match each other. 

For the outer expansion, diffusion can be ignored in equation (9). 
This leads to 


9n* 

9t’ 




t2 


(25) 


which is the well known equation for coagulation developed by Smolu- 
chowski and presented in Fuchs [1] . The solution for this equation is 


n’(t') 


1 

1 + yt’ 


(26) 


The inner expansion begins with the assumption that near the wall 
the only removal mechanism of any significance is diffusive deposition 
on the wall. It is further assumed that the particles are diffusing 
only along the x coordinate. These assumptions seem reasonable due to 
the extremely large* gradient which exists at the walls of the chamber. 
Thus the equation describing the depletion of the aerosol near the wall 
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is 


9n^ 

9t’ “ 9F^ 


( 27 ) 


where x is the distance from the wall. Carslaw and Jaeger [8] give a 
solution to this equation. 


n* (x’ ,t ’) « er£( 2 ^^. i ) 


( 28 ) 


This solution technique is presented as a possible alternative to 
the successive approximation technique presented earlier. However, 
the inner-outer expansion technique is a somewhat less desirable ap- 
proach. For the inner-outer expansion technique it is necessary to 
solve two problems while for the successive approximation technique 
only one problem need be solved. Further, matching the two solutions 
repeatedly is likely to be a problem whose difficulty is comparable to 
that of solving for the higher order approximations in the successive 
approximation technique. 
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CHAPTER III 


GRAVITY SOLUTION 


Equation 

A somewhat more difficult case is the solution of equation (12) in 
the presence of a gravitational field. A transformation is made which 
removes the sedimentation term. This transformation is due to Wilhelm 
[6] and is 

u(r',z',t') = f (r' ,z' ,t')exp(|z' - (29) 

This transformation yields 


9f 

3t' 


V^f - e exp( 




(30) 


The boundary conditions then become f(r',z',t') = 0 at the. wall and 
f(r',z',0) = exp(- l^')* 

Solution 

Equation (30) is solved by separation of variables in the same man- 
ner as was equation (15), resulting in 


f(r’.z'.t') = 


? Z A..J (X.r')sin(^l^^’)exp(-(X^ + -b ' ) ^ ) 

=1 i =1 ij o X n X n 


( 31 ) 


where 


/ ^ /^r ’ exp (- ) J (X , r * ) sin -a ' )dz ' dr ’ 

(X.r ’ )sin^ ( ^^^*z ')dz*dr* 

Q Q o 1 H 

4jTT(l - (-l)'^exp(- |l)) 

— n — 

X.Jj^CX.K^^ + j=^Tr^) 


Solution Near the Wall 

There are several areas in which this solution has difficulties. 

An examination of equation (29) will show where these difficulties oc- 

-3 

cur. Values for 3 range from 3.58 X 10 for small particles to 
12 

1.11 X 10 for large particles. Thus for certain particle sizes, the 
transformation given by equation (29) becomes so large that the results 
are meaningless. The solution does work for particles smaller than 
0.04 micrometers where sedimentation is much less important that dif- 
fusion, and for particles larger than 0.7 micrometers where diffusion 
is much less important than sedimentation and for long times when sed- 
imentation has essentially depleted the chamber. 

Results from the zero-gravity solution indicate that the effects of 
diffusive deposition are limited to a thin boundary layer. Thus dif- 
fusive deposition has a negligible effect on the concentration of the 
bulk of the chamber and it seems reasonable to expect that the problems 
of the above solution can be overcome by ignoring diffusion altogether. 
This leads to the equation 
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now applying the transformation in equation (11) the result is 


9u ^8u 

3t’ " 


(34) 


The solution to this is in terms of a unit step function, s(x). 
The unit step function is defined as follows 


s(x) = { 


0 

1 


X < 0 

X > 0 


So the solution to equation (34) is 


u(r',z',t') = s(|* - z' - gt') 


(35) 


This solution can be combined with the inner expansion solution of 
equation (28) to give the solution for the entire chamber. 
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CHAPTER IV 


RESULTS AND DISCUSSION 


Zero-Gravity 

Figure 1 gives the results of the settling velocity calculations 
for unit density spheres in air at 23C and 1 atmosphere. 

Figures 2-5 give the results of equation (23) which describes the 
depletion of a confined aerosol which is acted on only by coagulation 
and diffusive deposition. The first two figures give the concentration 
as a function of particle size for a location near the center of the 
chamber and times of 1 and 2 hours. These figures do not show any ef- 
fects of diffusive deposition on the aerosol concentration and are 
therefore identical to a curve for coagulation alone. Figures 4 and 5 
give similar information but for a location near the wall of the cham- 
ber. The dotted lines on the curves for the smaller particles show what 
the concentration would be if only coagulation were acting on the aero- 
sol at that point. Thus it can be seen that diffusive deposition is 
not very effective as a removal mechanism for the bulk of the chamber. 
Near the wall, however, the situation is reversed. There is a steep 
gradient in concentration which drives diffusive deposition. Thus for 
a thin boundary layer the dominant removal mechanism is diffusive de- 
position. The steep gradient near the wall will reduce equation (14) 
to a one dimensional equation except at the corners. The solution to 
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Figure 1. PARTICLE SETTLING VELOCITY. 





Figure 3. PARTICLE DEPLETION IN THE ABSENCE OF GRAVITY: IMPORTANCE OF COAGULATION. 



p IP 



Figure 5. PARTICLE DEPLETION IN THE ABSENCE OF GRAVITY: IMPORTANCE OF COAGULATION. 
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the one dimensional diffusion equation (28). Figures 6 and 7 are con- 
centration profiles for the region near the wall of the chamber and 
show the extent of the action of diffusive deposition for times of 1 
and 2 hours respectively. Thus for the zeroth approximation, equation 
(15) gives the concentration for the entire chamber. 

The zeroth approximation is sufficiently accurate for all but the 
most demanding applications. The higher order solutions are very dif- 
ficult to obtain and contain only a small correction to the solution 
given by equation (23). Table I gives the values of the terms of equa- 
tion (14) for various locations within the chamber and times of 1 and 2 
hours. The values for locations near the center of the chamber are 
large relative to the other terms in equation (14) but are quite small. 
The size of e for these locations is a further indication that there 
are no gradients and thus diffusive deposition is not important. The 
values are so large because of truncation error in evaluating equation 
(23) and not any fluctuations in aerosol concentration at those points. 
If the truncation error were not present it is expected that the value 
of e would be zero. Substituting e = 0 into equation (14) for the first 
approximation would yield exactly the same result as the zeroth approx- 
imation. Thus for those situations where e = 0 the zeroth approxiamtion 
is exact. 

The inner-outer expansion solution embodied in equations (25) 
through (28) give further evidence that the zeroth approximation is ac- 
curate because this solution technique gives the same results as those 
of equation (15). That one dimensional diffusion is the dominant re- 
moval mechanism near the wall is also shown by the fact that the solu- 
tion to equation (27) predicted the difference between the curves for 







o 0.4 



Unit density shperes in still air at 23C and 1 atm. 
Cylindrical chamber: radius = 25cm. height = 50cm. 
Location; z “ 25cm. t = 2hr. 
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Table I 


VALUES OF TERMS IN EQUATION (14) 


R(cm) 

r (cm) 

t(hr) 

3u 

3t' 

V^u 

e 


lE-6 

5 

1 

0. 

3.05E-6 

3.05E-6 

lE-6 

22.5 

1 

18.96 

18.97 

5.80E-3 

lE-6 

5 

2 

0. 

4.18E-16 

4.18E-16 

lE-6 

22.5 

2 

87.69 

88.03 

0.340 

lE-5 

5 

1 

0. 

4.09E-9 

4.09E-9 

lE-5 

22.5 

1 

0. 

4.09E-9 

4.09E-9 

lE-5 

5 

2 

0. 

2.26E-12 

2.26E-12 

lE-5 

22.5 

2 

0. 

2.26E-12 

2.26E-12 

lE-4 

5 

1 

0. 

5.18E-6 

5.18E-6 

lE-4 

22.5 

1 

0. 

5.40E-6 

5.40E-6 

lE-4 

5 

2 

0. 

6.47E-6 

6.47E-6 

lE-4 

22.5 

2 

0. 

6.48E-6 

6.48E-6 

lE-3 

5 

1 

0. 

7.46E-6 

7.46E-6 

lE-3 

22.5 

1 

0. 

1.35E-5 

1.35E-5 

lE-3 

5 

2 

0. 

1.41E-5 

1.41E-5 

lE-3 

22.5 

2 

0. 

2.24E-5 

2.24E-5 
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coagulation alone, depicted by the dotted lines, and the curves for 
both coagulation and diffusive deposition in figures 4 and 5 to within 
1 %. 

Table I shows that the values for e for locations near the wall are 
quite large relative to those far from the wall. Comparison with the 
values of the other terms in equation (14) shows that they are negli- 
gible. Thus the approximation calculated by ignoring e for those re- 
gions was good. 

Gravity 

Equation (33) is the model of a confined aerosol which is acted on 
by sedimentation as well as coagulation and diffusive deposition. Fig- 
ures 8 and 9 give the results of this model. Figures 8 and 9 are 
graphs of the concentration as a function of particle size for a loca- 
tion near the bottom of the chamber. An inportant result of this model 
is the sharp gradient which appears in the vertical extent of the cham- 
ber as a result of the action of gravity. This gradient moves downward 
with a velocity equal to the terminal settling velocity of the particles. 
Above the line of the gradient there are no particles. Below the line 
the concentration changes as if the particles were being acted on by 
coagulation and diffusive deposition. 



Figure 8. PARTICLE DEPLETION IN THE PRESENCE OF GRAVITY AND 


Unit density spheres in still 
air at 23C and 1 atm. 

Cylindrical chamber : radius 

= 25cm. height = 50cm. 

Location: R = 5cm. z = 5cm. 

t = Ihr. 


NJ 


■6 -t. I « I ^ 
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COAGULATION. 




Figure 9. PARTICLE DEPLETION IN THE PRESENCE OF 



^VITY AND COAGULATION. 
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CHAPTER V 


CONCLUSIONS 


The aerosol depletion model given by equation (14) represents the 
action of a confined aerosol both in the presence and absence of gravity 
for the entire chamber. It was found that in the presence of a gravi- 
tational field, sedimentation was the dominant mechanism, even for small 
particles. When there is no gravitational field, coagulation is the 
most significant mechanism for removing particles from the bulk of the 
chamber . 

An advantage of the zero-gravity model is that it embodies the so- 
lution for the entire chamber. The more conventional technique of in- 
ner-outer expansions requires that two different problems be solved and 
the solutions then matched. The matching problem is likely to present 
difficulties comparable to those of solving for the higher order approx- 
imations. 

The results presented are those of the zeroth approximation. 

However, the higher order approximations would only yield a small cor- 
rection to the zeroth approximation and thus may not be worthwhile in 
most cases. 

The situation in the presence of gravity yields an important re- 
sult in the sharp boundary which moves downward through the chamber. 

This result is not unreasonable and is substantiated by the results of 
the one dimensional diffusion work. This work indicated that aerosol 
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particles do not diffuse very far from the location of the gradient 
which drives the diffusion. Nonetheless, diffusion is present and the 
boundary would be more or less blurred by it . 

The model is not complete in that it does not directly deal with 
the shifting size distribution of particles caused by coagulation. 

Fuchs [1] indicates that there has been no significant time dependence 
in the coagulation coefficient detected in the terrestrial experiments 
done to date. This is because the expected increase in the coagulation 
coefficient due to the increased particle size is almost exactly offset 
by the decrease in diffusivity of the larger particle. The diffusion 
coefficient for the new particle would be less, but this would have the 
effect of decreasing the effect of diffusive deposition on the particle 
concentration. Thus this effect should not require a major correction. 
The terminal settling velocity of the new particles would be greater 
and this could cause a problem for the gravity model. It is expected 
that for the lower initial concentrations and short times that the re- 
sults of the gravity model would be valid because there would only be a 
few of these larger particles formed. An attempt to deal with this 
problem would follow the lines of research presented by Tolfo [9]. This 
would Involve writing a model such as this for each size class and sol- 
ving them simultaneously. Some further work that would be instructive 
would be to find the higher order approximations, although it is ex- 
pected that they will not provide much additional information. 

There is much experimental work that needs to be done. Data are 
needed to validate this model. These data should take the form of 
aerosol concentrations for various times and locations within the cham- 
ber, both in the presence and absence of gravity. The values of the 
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coagulation coefficient and diffusion coefficient used in this model 
were theoretical. A zero-gravity experiment would be instrumental in 
obtaining measurements of these quantities* These measurements are 
severely hampered in terrestrial experiments by the action of gravity. 
Zero-gravity experiments could be designed which would provide direct 
measurements of these quantities, thus validating that portion of this 
mode^ and advancing our knowledge of aerosol science. 
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APPENDIX A 
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PROGRAM SOLUTN ( INPUT f OUTPUT , I aPES= INPUT t TAPEb=OUT PUT ) 

SOL 

SOL 

10 

20 

• • 


• SOL 

30 


PROGRAM CALCULATES THE ZEROTH APPROXIMATION TO THE SOLUTION 

SOL 

40 


OF The AEROSOL DECAY EQUATION FOR A ZERO GRAVITY ENVIRONMENT 

SOL 

50 

• • 

TN'THIS PART variables ARE INITIALIZED AND PROGRAM OPERATION 

» SOL 

50 


SOL 

70 


IS CONTROLLED 

SOL 

80 

» • 


* SOL 
SOL 

90 

100 


dimension f ( 11 * 1 1 ) .u( 1 1 1 n ) »epslon( ii • i n *z < n i «concn( iif in • 

SOL 

no 

] 

1 SEBSlN n 1 > «nUDZ ( 1 1 • U ) tOUOH (Util) f SERJO (11) 

SOL 

120 


.COMMON/ZZ/BETA.PI ,PI02.PI02SQfBTASQtEXPBTA 

SOL 

130 


COMMON/RH/TIHE 

SOL 

SOL 

140 

150 

* • 

* * * * INITIALIZE 

SOL 

SOL 

160 

170 


TlMPRM = 3f>00. 

SOL 

180 


CmMHA0=25. 

SOL 

190 


CHKwT=50. 

SOL 

200 


CONCl=*i ,E4 

SOL 

210 


ip=n 

SOL 

220 


IZ=1 1 

SOL 

230 


RINC=. 1 

SOL 

240 


Z1NC=.Z 

SOL 

250 


PI=4.®ATAN(1.) 

SOL 

260 


PI02=PI/2. 

SOL 

270 


PI02SQ=P102*PI02 

SOL 

280 


ITEWM2 = 3rio 

SOL 

290 


IT£RmR=500 

SOL 

300 

10 

CONTINUE 

SOL 

310 


IPASS=«l 

SOL 

320 


READ <5f 20) PfVfCOAG.D 

SOL 

330 

?0 

FORMAT (4F10.0) 

SOL 

340 


IF tEOF (S) ) 270f 30»270 

SOL 

350 

30 

CONTINUE 

SOL 

360 


WRITE (6*40) RfVfDiCOAGtTlMPRM.CONCI 

SOL 

370 

40 

FORMAT (lHlf//*10Xf IBHTHIS OUTPUT IS FOR#//flbXt 13HPARTICLE 5I2E50L 

380 


1 . 15X, IPEe . 1 I /| 15X* 17H5ETTLING VE LOC I T Y f 1 1 X . 1 PE H . A . / i 1 bX t tllHOIFFSOL 390 
2USI0N COEFF ICIENT ,7X . IPE n. A,/, 15X. 23HC0AGUL A r 1 ON COE FF I C I ENT . S X . SOL AOO 



31PEll.<t,/,15Xi AHTIMEt2AX,0PF8.U/«15Xt 

21HINITIAL CONCENTRAT lON.SOL 

410 


47X» |PC8. 1 f //) 

SOL 

420 


0ETA=CHMPAD«V/D 

SOL 

430 


GAMHA=CHMRA0®CHMRAD»C0AG*C0NCI /D 

SOL 

440 


T IME = T 1MPRM«0/ (CHHRAO®CHMRAD) 

SOL 

450 


BTAS0 = BETA*Bfc‘.TA/4. 

SOL 

460 


BTAOV2=RETA/2. 

SOL 

470 


F.XPRTA = EXP (-PETA^CHMHT / (2,*CHMRAO) ) 

SOL 

480 


BTaTIM=»TA50«TIME 

SOL 

490 


XPTIME=EXP(BTATIM) 

SOL 

500 


ghatih=gamma«time 

SOL 

510 


T0GMAT=2,*GMATIM 

SOL 

520 

C 


SOL 

530 

C 


SOL 

540 

C 


SOL 

550 


Z ( 1 ) = 0 . 

SOL 

560 


DO 50 I=2fIZ 

SOL 

570 


Z ( I) =Z ( I-l ) ♦ZINC 

SOL 

580 


50 CONTINUE 

SOL 

590 


00 60 1 = 1 *IZ 

SOL 

600 


CALL SINES (Z (I ) tSEMSIN(l) tlTERMZ) 

SOL 

610 


60 CONTINUE 

SOL 

620 

c 


C r.i 

6ZD 

c 


SOL 

640 

c 


SOL 

650 


RAD:^0. 

SOL 

660 


DO 70 I=l*I« 

SOL 

670 


CALL BESSEL ( RAO • 5EH JO ( I ) t I TERMH ) 

SOL 

680 


RAO=RAO«RINC 

SOL 

690 


70 CONTINUE 

SOL 

700 

c 


SOL 

710 

c 


SOL 

720 



non orjr> ooo 
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SOL 

730 


00 no 1>1»IH 

SOL 

740 


00 100 J=l»I2 

SOL 

750 


F (I t J) =StHj0 ( I ) *SEPSIN( J) 

SOL 

760 


XEXP=BTA0V2*/ ( J)-BTATIM 

SOL 

770 


IF (XFXP,GT.600.) GO TO 80 

SOL 

7B0 


W(ltJ)=F(I»J)*EXP(XEXP) 

SOL 

790 


CONCN<I*J)=U(I.J)/(l.»eMATIM«U(ItJ)) 

SOL 

800 


GO TO 100 

SOL 

810 

80 

CONTINUE 

SOL 

820 


WRITE <6f90) I f JfF ( I f J) tXEXP 

SOL 

830 

90 

format <l5Xf AHI a fl3»3X* AHJ » fI3*3X» 9HFtI#J) 

* ,1PE11.4,3XS0L 

840 


It 7HXEXP = tlPEll.4) 

SOL 

B50 


IPASS=-I 

SOL 

860 

100 

continue 

SOL 

870 

110 

CONTINUE 

SOL 

BBO 


IF UPASS.LT.O) go to 200 

SOL 

890 



SOL 

900 

# • 

* • • • CALCULATE DEHIVATIVES FOH EPSILON 

SOL 

91 0 



SOL 

920 


DO 120 I = I , IR 

SOL 

930 


OUDR (1 , I ) = (U( 1 • I )-U(2t 1 > )/ (-RING) 

SOL 

940 

120 

CONTINUE 

SOL 

950 


DO 149 1=2, IR 

SOL 

960 


DO 130 J=mz 

SOL 

970 


OUDR 1 I , J) = (U( 1-1 , J)-U(I t J) ) /(-RING) 

SOL 

980 

130 

CONTINUE 

SOL 

990 

1 AO 

CONTINUE 

SOL 

1000 


DO 150 1 = 1, IZ 

SOL 

1010 


DUD2 (I • 1 )= (U(I » 1 ) -U( I f2) ) / <-ZINC) 

SOL 

1020 

150 

CONTINUE 

SOL 

1030 


00 1 70 1 = 1 , IR 

SOL 

lOAO 


DO 160 J=2tlZ 

SOL 

1050 


OUDZ(I*J)s(U(IfJ-l)-U(ItJ) )/(-ZINC) 

SOL 

1060 

160 

CONTINUE 

SOL 

1070 

170 

CONTINUE 

SOL 

1080 



SOL 

1090 

« • 

* * * * CALCULATE EPSILON 

SOL 

1100 



SOL 

1110 


DO 190 1=1, IR 

SOL 

1120 


DO IPO J=1,IZ 

SOL 

1130 


EPSLON ( I . J) = ( TOGMAT/ ( 1 .«GMaTIM*U ( I , J) ) ) * (DUOR ( I 

f J) *DUOR ( I • J) SOL 

1 lAO 

1 ♦OnoZUf J)«0U0Z(I#J) ) 

SOL 

1150 

160 

continue 

SOL 

1160 

190 

CONTINUE 

SOL 

1170 



SOL 

1160 

• • 

* • » • OUTPUT variables 

SOL 

1 190 



SOL 

1200 

200 

continue 

SOL 

1210 


►RITE (6,210) 

SOL 

1220 

210 

format ( 15X* 5HCONCN) 

SOL 

1230 


call terms ( IR, IZtCONCN) 

SOL 

1240 


WRITE (6f220) 

SOL 

1250 

220 

format ( isx » 7HEPSILON) 

SOL 

1260 


CALL terms (!«♦ I/f EPSLON) 

SOL 

1270 


WRITE (6#230) 

SOL 

1280 

230 

FORMAT ( 15X# IMU) 

SOL 

1290 


CALL TERMS (IR»I2*U) 

SOL 

1300 


write (6,240) 

SOL 

1310 

240 

FORMAT (ISXt 4H0UDR) 

SOL 

1320 


call terms (IR,IZ,DUDR) 

SOL 

1330 


WRITE (6,250) 

SOL 

1340 

250 

format nsx. AHDUDZ) 

SOL 

1350 


call terms (IR,IZ,DU0Z) 

SOL 

1360 


WRITE (6,260) 

SOL 

1370 

260 

format nsxt IHF) 

SOL 

1380 


call terms (IR,IZ,F) 

SOL 

1390 


GO TO 10 

SOL 

1400 

270 

CONTINUE 

SOL 

1410 


STOP 

SOL 

1420 



SOL 

1430 


END 

SOL 

1440 
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SUBROUTINE BESSL (GNUtXtBESSEU BES 

BtS 
hES 

THIS 5R CALCULATES THE BESSEL FUNCTlONf J(X)t OF ORDER GNU BES 

«•••«•»«•««««»* «*•»•»*»*«•«• hES 

BES 

DATA AA21/.25/iAAA 1/-.375/»AA42/. t5b25/»AA6l/l.B7b/f BES 

1 AA62/-1 .lS62S/.AA63/.U7lB75/,AA8l/-l9.6875/tAAB2/14.2265625/f BES 

2 AAB3/-2.3867lB75/f AA84/,09b2l48437/*AAl0X/3S4.3Tb/t BES 

3 AAl02/-'277.8Tb/f AAl03/bB«2246093h/f AAIOA/-^. i00bttS938/* BES 

4 AAl0a/.0809326l7l/#BB2l/.b/*BB4l/-.25/»BH42/.U4l666b66666/f BES 

5 P061/,75/*HBb2/-.35/»BB63/.012S/.»HH8l/-5.62b/fBB82/3«6026TB8571/«HES 

6 8B83/-.424107I42q/,BB84/b.b803S7l43/*BBl01/78. ?b/ t BB 1 02/-b8 • / t BES 

7 BD I 03/ 10, 28645833/ #BB I 04/-. 347222222222/ »BB I Ob /3, 0381 9444 4E-3/ RES 

HES 

• ♦ • * INITIALIZE HES 

HES 

PI=4.*ATAN(1.) BES 

P102=Pl/2. BES 

ALPHA=GNU*6NU-. 25 BES 

T=l,/X BES 

TSQoT^T BES 

RES 

^ • CALCULATE COEFFICIENTS BES 

BES 

A2=ALPHA»AA21 BES 

A4=tAA42«ALPHA*AA41) • ALPHA BES 

A6 = ( ( AA6 3*^ALPHA-»AA6 2) »ALPHA*AA61 > •ALPHA BES 

AB=f{(AA840ALPHA^AA8 3)»ALPHA4AA82) *»ALPHA4AA8l ) * ALPHA BES 

AlO=t H ( (AAl 05* ALPHA* AA104) *ALPHA*AAl 03) •ALPHAS AAIU2) •ALPHA^AAlOl ) •BES 
lAl.PMA HES 

B2=BR21*ALPHA BES 

fl4=1BR42«ALPHA*H04l ) "ALPHA BES 

B6= f (HHb 3 "ALPHA *8862) "ALPHA^BBbl ) ♦ ALPHA BF.S 

R8 = f ( (Bnn4"ALPHA*BBB3) **ALPHA*BB82) "ALPHA *0081 ) "ALPHA BES 

01 03 ( ( ( (BR10 5"ALPHA*BB104 ) " ALPH A ♦ HB 1 0 3 ) * ALPHA * BB i 02 ) "ALPHA*B0 1 0 1 ) "BES 
lALPHA BES 

BES 

► * • * * CALCULATE INTERMEDIATE FUNCTIONS BES 

BES 

B3(<((A10"TSQ*A8)*T5U^A6)*TSQ*A4)»TSQ*A2)*TSO*1« HES 

PRTPHl=((((B10"TSQ*B8)*TSQ*U6)*TSQ*B4)*TSU*H2)*TSO*l. BES 

PHI3PRTPHI/T- (GNU^.5) "PI02 BES 

BES 

I • • » • CALCULATE BESSEL FUNCTION HES 

HES 

BESSEL3B*S0RT (T/PI02) "COS (PHI ) BES 

RETURN BES 

C HES 

END HES 


10 

20 

30 

40 

50 
60 
70 
60 
90 
100 
no 
120 
130 
140 
ISO 
160 
170 
160 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
42 0 
430 
440 
450 
460 
470 
480 
490 
500 
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SUBROUTINE BLAMOA ( XL AMOA » I TERM ) 


FINDS The first ITERM ROOTS OF THE 2ERO ORDER BESSEL FTN 


"Dimension xlamda(iterm) 

initialise 

PI = <».*ATANn.) 

A = 2.4 

ERWOR=l ,F-8 
DO 60 I=l*ITERM 
B = A^l . 

X = B 

IF (A.LE.IO.) call HESJ I A * 0 • X JOF A t U E-6 • I EH ) 
IF (B.LE.IO.) CALL BESJ ( P • 0 . X JOF0 * 1 .E-6 * lER ) 
IF (A.DT.10.) call BESSL (0.fA-»XJOFA) 

IF (B.GT.IO,) CALL BESSL (O.tB^XJOFB) 

10 continue 

IF (X.GT.IO.) CALL BESSL (O.tX^XJOFX) 

IF (X.LE.IO.) CALL BESJ ( X t 0 • XUOF X # I , E -6 f I EH ) 
IF (XJOFA*XJOFX) 20ib0i30 

20 continue 

QrX 

XJ0FB=XJ0FX 
GO TO 40 
30 CONTINUE 

A = X 

X J0FA=XJ0FX 
40 continue 

IF (ARS(XJOFX) ,LE. ERROR) GO TO 50 
X= ( A«XJOFB-B<>XJOFA) / ( X JOFB-X JOF A ) 

GO TO 10 
SO CONTINUE 

A = XO, 

XLAMOA u ) =X 
60 CONTINUE 
RETURN 

END 


BLA 10 
BLA 20 
HLA 30 
RLA 40 
HLA so 
BLA 60 
BLA 70 
BLA 80 
BLA 90 
BLA 100 

RLA no 

BLA 120 
BLA 130 
BLA 140 
BLA ISO 
BLA 160 
BLA 170 
BLA 180 
BLA 190 
RLA 200 
BLA 210 
RLA 220 
BLA 230 
BLA 240 
BLA 250 
HLA 260 
HLA 270 
BLA 260 
HLA 290 
BLA 300 
BLA 310 
RLA 320 
HLA 330 
BLA 340 
HLA 350 
BLA 360 
HLA 370 
BLA 380 
BLA 390 
BLA 400 
BLA 410 
BLA 420 


SUBROUTINE SINES IZiFilTERM) 


THIS SR CALCULATES THE SIN SERIESt F* FOR POSITION Z 


DIMENSION A (350 ) »SINE(350 ) t TERM (350 ) 

COMNON/RR/T IME 

COHMON/ZZ/BETA»PI #PIO2tPl02SQ*BTA9Q»EXPBTA 

F = o: 

XMNl=l, 

00 10 ITRRrlflTERM 
XMNl = XMNl (-1 . ) 

XTWM=FL0AT (ITRM) 

tXPTlM = r XP (-XTRM*XTHM*PI02SCJ«TIME) 

AR0sXTRH«PI02*Z 

A ( ITRM) SXTHM*P 102* ( 1 • -XMN 1 *EXPBT A ) / ( BT ASQ ♦ X TRM*XTRM*P I025Q ) 
SINE(ITRM)=SlN(AHG) 

TERM(nRM)=A(nHM)*SlNE(ITHM)<»EXPTIM 
F=F*TERM( iTRM) 

10 CONTINUE 
RETURN 
C 

END 


SIN 10 
SIN 20 
SIN 30 
SIN 40 
SIN 50 
SIN 60 
SIN 70 
SIN 80 
SIN 90 
SIN 100 
SIN 110 
SIN 120 
SIN 130 
SIN 140 
SIN 150 
SIN 160 
SIN 170 
SIN 180 
SIN 190 
SIN 200 
SIN 210 
SIN 220 
SIN 230 
SIN 240 


o o o o n 
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SUBROUTINE BESSEL (R#Ff ITEPM) BSS 

HSS 

ess 

THIS SH CALCULATES THE JU StRIESf Ft FOR POSITION R BSS 

BSS 

BSS 

-^lHt>jSTON XLHDA(500 ) tBF.SSLI (SOO ) tBESSLO (SOU I tTERHISOO ) t BSS 
1 A(S00 ) HSS 

COMMON/RR/TIME HSS 

F=0i BSS 

ERROR^l . e-6 HSS 

CALL BLAMDA ( XL MDA 1 1 TERM ) BSS 

DO 80 1TRM=I . ITERM BSS 

IF (XLMOAdTRM) •GT.IO.) GO TO 20 BSS 

CALL BESJ (XLHDA(ITRM) # 1 f BES5L I ( I TRM ) f ERROR* I EH) BSS 

IF (lt:R.NE«0) WRITE (btIO) 1 EP « XLHDA ( 1 TRM ) iHESSLl (I TRM) BSS 

10 format (10X» 5HIER= * 13 % 9H LAMOA a ♦IPEll.At bH Jl » *IPE11.4)BSS 
GO TO 30 BSS 

20 CONTIMUF BSS 

CALL HESSL < I • • XLMD A ( I T RM ) * BESSL l( I TRM ) ) HSS 

30 CONTINUE BSS 

ARO = R»XLMr)A ( ITRM) HSS 

IF (APG.EO.O,) GO TO 60 HSS 

IF (ARG.GT.10.) GO TO bO P5S 

CALL PESJ ( ARGf 0 *BESSL0 ( I TRM) *ERROR. lER) BSS 

IF (lER.NE.O) WRITE ( 6 * AO) I ER * ARG * BE SSL 0 ( I TRM ) BSS 

40 format noX. 6HIER = *13* 7H ARG = *lPEU.4i 6M JO = tlPEll.4) BSS 
GO TO 70 BSS 

50 continue BSS 

CALL PESSL (0. *AR6fBESSL0 ( ITRM) ) RSS 

GO TO 70 BSS 

60 CONTINUE BSS 

BESSLO (ITRM) si • BSS 

70 CONTINUE BSS 

BO CONTINUE BSS 

00 90 1=1 *1TERM BSS 

A ( I ) =2./(BtSSLl (I>*XLMDA(1) ) BSS 

90 CONTINUE HSS 

DO 100 ITRM=1*ITERM BSS 

CXPTIH=EXP(-XLMDA( ITRM) *XLMDA ( I TRM) *T IME) BSS 

TFHM(ITRH)=A(I TRM) *BES5lO ( ITRM) «EXPTIM BSS 

F=FdFRH(lTRM) BSS 

100 CONTINUE BSS 

return RSS 

BSS 
BSS 


10 

20 

30 

40 

SO 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

34() 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 


C 


END 
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SUBROUTINE TERMS ( I TERM . JTEHM t A ) 
DIMENSION A ( ITERM»JTERM) 

WRITE (6*10) 

10 format (/) 

DO 30 K=1*JTEHM 

W45ITE- (6*20) (A(J»K) f 

?0 format (3Xf n (1PE9.2»3X) ) 

30 CONTINUE 

IF flTERM.LE.ll) GO TO 70 
WRITE (6f 10) 

DO ^0 K=1 t JTERM 

WRITE <6»20) (A(J»K)tJ=12»22) 

40 CONTINUE 

IF « lTtRM.LE.22) GO TO 70 
WRITE <6. 10) 

DO 50 K=1*JTEWM 

WRITE (G»20) (A(J»K)tJ=*2J,33) 

50 continue 

IF UTERM.LE.33) 60 TO 70 
WRITE (6f 10) 

DO 60 K=1 . JTERM 

WRITE (6»20) (A ( JtK) • J=34t44) 
60 CONTINUE 
70 CONTINUE 

WRITE (6f0O) 
flo format (//) 

RETURN 

END 


TER 10 
TER 20 
TEH 30 
TER 40 
TER 50 
TER 60 
TER 70 
TEH 80 
TEH 90 
TER 100 
TER 110 
TER 120 
TEH 130 
TER 140 
TEH 160 
TER 160 
TEH 170 
TEH 100 
TEH 190 
TER 200 
TEH 210 
TEH 220 
TER 230 
TER 240 
TER 250 
TER 260 
TER 270 
TER 280 
TEH 290 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM SOLUTN ( INPUT tOUTPUT . TAPE5= INPUT t TAPE6*0UTPUT ) 


PROGRAM CALCULATES THE ZEROTH APPROXIMATION TO THE SOLUTION 
OF THE AEROSOL DECAY EQUATION FOR A GRAVITY ENVIRONMENT 

IN THIS PART VARIABLES ARE INITIALIZED AND PROGRAM OPERATION 
IS CONTROLLED 


DIMENSION F < n « 1 1 ) tU< 11 « 1 1) tEPSLONU 1 *1 1 ) t Z < 1 1 ) tCONCNIIl 1 1 i) » 

I SERSIN ( 1 1 ) «DUOZ (11*11) *DUOR (11*11) tSERJO (11) 
COMHON/ZZ/nETAfPI*P102fPI02SQ*BTASQ*EXPBTA 
COMMON/RR/TIHE 

^ • « • « » INITIALIZE 

T IMPRM=7200. 

COMC I»1 ,E5 
CHMR4D»25. 

CHMHTobO. 

H=CHM^T/CHMRAD 
IP = 1 I 
IZ = 1 1 
RINC=.l 
ZINC=.2 

PI=A.»ATAN( 1 , ) 

Pl02=PI/2. 

PI02SQ=PI02*PI02 

ITERH2=350 

ITe:Rm»s500 

Z(1)=0, 

DO 10 l^ZflZ 

Z(I)=Z(I-1) ♦ZINC 
10 CONTINUE 
20 CONTINUE 

READ (5*30) Rf V*C0AG*0 
30 FORMAT (4F10.0) 

IF (EOF (S) ) 210*40*210 
40 CONTINUE 

WRITE (6*50) W* V *0*C0AG*T IMPRM,C0NCI 
50 FORMAT (iHlt/ZtlOXt IflHTHIS OUTPUT IS F0R*//*l5Xi 13HPARTICLE SIZ( 
l*15X*iPFa.l*/*15X* 17HSETTLING VELOC I T Y * 11 X . 1 PE 11 . 4 , / * 1 bX . 21H0IFI 
2USI0M COEFFICIENT *7X. IPEll .4,/, IbX* 23HC0AGUL A T I ON COEF F I C I ENT • bX < 
31PEU.4*/f 15X* 4HTlME*24XfOPF0,lfc/f 15X. ZIHINITIAL CONCENTRATION! 
47X* IPEB. I •//) 

RET A=QHMRAD*V/D 

gamma=chmmao*chmrad*coag*conci/d 

TIMEsT IMPRM«0/ (CHMRAD«CHMRAD) 

BTAS0iiHETA*HETA/4. 

HTA0V2rRf:TA/?, 

EXPOTA=rxP (-HETA^CHMHT/ (2,*CHMRAO) ) 

btatim=rtasq*time 

HTTlMF.rHfTA<*TlHE 

XPTlME=rxP(BTATIH) 

GMAT IN = GAMMA«»T IME 
T0GMAT=2,*GMATIM 

••••«* CALCULATE TOTAL SOLUTION 

no 7 0 I -1 » IP 
00 60 J=l*IZ 

STPFNI=Z ( J) -BTT IME 

STPFN2=STPFN1-H 

IF (ST PEN 1 .LT *0.) SPFNI«0, 

IF (STPFNl .GE.O. ) SPFNi=l. 

IF (STPFN2.LT.0. ) SPFN2sO. 

IF (STPFN2.GE.0. ) SPFNZf 1 • 
un • J)=SPFN1-SPFN2 

CONCN(I*J)sU(I* J)/(1.^GMATIM*U(I .J> ) 

60 continue 

70 CONTINUE 


SOL 

10 

SOL 

20 

SOL 

30 

SOL 

40 

SOL 

50 

SOL 

60 

SOL 

70 

SOL 

60 

SOL 

90 

SOL 

100 

SOL 

110 

SOL 

120 

SOL 

130 

SOL 

UO 

SOL 

150 

SOL 

160 

SOL 

170 

SOL 

180 

SOL 

190 

SOL 

200 

SOL 

210 

SOL 

220 

SOL 

230 

SOL 

240 

SOL 

250 

SOL 

260 

SOL 

270 

SOL 

260 

SOL 

290 

SOL 

300 

SOL 

310 

SOL 

320 

SOL 

330 

SOL 

340 

SOL 

350 

SOL 

360 

SOL 

370 

SOL 

380 

SOL 

390 

SOL 

400 

SOL 

410 

;SOL 

420 

SOL 

430 

SOL 

440 

SOL 

450 

SOL 

460 

SOL 

470 

SOL 

460 

SOL 

490 

SOL 

500 

SOL 

510 

SOL 

520 

SOL 

530 

SOL 

540 

SOL 

550 

SOL 

560 

SOL 

570 

SOL 

580 

SOL 

590 

SOL 

600 

SOL 

610 

SOL 

620 

SOL 

630 

SOL 

640 

SOL 

650 

SOL 

660 

SOL 

670 

SOL 

66() 

SOL 

690 

SOL 

700 

SOL 

710 

SOL 

720 
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DO ao I»1.U SOL 

U (I « 1 1 ) kO . SOL 

unit n-o. SOL 

CONCNUtin^O. SOL 

coNCNni*n»o, sol 

80 CONTINUE SOL 

SOL 

CALCULATE DERIVATIVES FOR EPSILON SOL 

SOL 

DO 90 1=1 »IR SOL 

DUORdf I)*(UUtI)“U(2f I) )/(-RINC) SOL 

90 CONTINUE SOL 

00 no 1 = 2. IP SOL 

DO 100 J=1#I2 SOL 

DUOR (I.J)=IU(I-1.J>-U(I»J))^(-RINC) SOL 

100 continue sol 

no CONTINUE SOL 

DO 120 I=l • IZ SOL 

DUDZ { I f 1 )= (U(I . 1 )-U( I .2) )/ (-ZINC) SOL 

120 CONTINUE SOL 

DO UO I = l,IR SOL 

DO 130 J=2.IZ SOL 

DUDZ ( I ♦ J) = (U( I • J-1 )-U( I f J) ) Z (-ZINC) SOL 

130 CONTINUE SOL 

UO CONTINUE SOL 

SOL 

CALCULATE EPSILON SOL 

SOL 

00 160 I = niR SOL 

DO ISO J=1 *IZ SOL 

IF (ABS (DUDR (I t J) ) • GE • 1 • E27 0 . OR . ABS ( OUOZ (1 1 J) ) .GE. 1 .£270 ) G050L 
I TO ISO SOL 

EPSLON ( I . J) = (TOGMAT/ ( 1 . ♦GMAT IM»U (I . J) ) ) <* (DUDR ( I • J) *DU0R ( I . J) SOL 
1 ♦OUnZUf J)*OUOZ(I#J) ) SOL 

150 CONTINUE SOL 

160 CONTINUE SOL 

SOL 

••«••• OUTPUT variables SOL 

SOL 

CALL terms ( IR. IZ.CONCN) SOL 

WRITE (6# 170) SOL 

170 format (ISX* 7HEP5ILON) SOL 

CALL terms (IR.lZfEPSLON) SOL 

WRITE (6. 180) SOL 

180 format nsx. iHU) sol 

CALL terms (IRfIZfU) SOL 

WRITE (6.190) SOL 

190 format USX. 4HDUDR) SOL 

call terms (IR.IZ.OUDR) SOL 

WRITE (6.200) SOL 

200 format (ISX. 4H0UDZ) SOL 

CALL TERMS (IR.IZ.OUDZ) SOL 

GO TO 20 SOL 

210 CONTINUE SOL 

STOP SOL 

C SOL 

FNH sni 


730 
740 
750 
760 
710 
780 
790 
BOO 
810 
B20 
B30 
64 0 
B50 
B60 
870 
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B90 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1000 
1090 

noo 
1110 
1120 
1130 
1 140 
1 150 
J 160 
1170 
1160 
1190 
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1230 
1240 
1250 
1260 
1270 
1260 
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